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Abstract. We have measured the effective mass, m, and Lande g factor in very 
dilute two-dimensional electron systems in silicon. Two independent methods 
have been used: (i) measurements of the magnetic field required to fully polarize 
the electrons' spins and (ii) analysis of the Shubnikov-de Haas oscillations. We 
have observed a sharp increase of the effective mass with decreasing electron 
density while the g factor remains nearly constant and close to its value in bulk 
silicon. The corresponding strong rise of the spin susceptibility x <^ 9™ ma^y 
be a precursor of a spontaneous spin polarization; unlike in the Stoner scenario, 
it originates from the enhancement of the effective mass rather than from the 
increase of g factor. Furthermore, using tilted magnetic fields, we have found 
that the enhanced effective mass is independent of the degree of spin polarization 
and, therefore, its increase is not related to spin exchange effects, in contradiction 
with existing theories. Our results show that the dilute 2D electron system in 
silicon behaves well beyond a weakly interacting Fermi liquid. 
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1. Introduction 

At sufficiently low electron densities, two-dimensional (2D) electron systems become 
strongly correlated, because the kinetic energy is overpowered by energy of electron- 
electron interactions. The strength of the interactions is usually characterized by 
the ratio between the Coulomb energy and the Fermi energy, = E^/Ep, which, 
assuming that the effective electron mass is equal to the band mass, in the systems 
with single- valley spectrum reduces to the Wigner-Seitz radius, l/lirnsY^^aB (here rig 
is the electron density and as is the Bohr radius in semiconductor). There are several 
suggested candidates for the ground state of the system, for example, (i) Wigner 
crystal characterized by spatial and spin ordering (ii) ferromagnetic Fermi liquid 
with spontaneous spin ordering [2], and (iii) paramagnetic Fermi liquid [3]. In the 
strongly-interacting limit (r^ ^ 1), no analytical theory has been developed to date. 

* To whom correspondence should be addressed 
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According to numeric simulations 0], Wigner crystallization is expected in a very 
dilute regime, when reaches approximately 35. The refined numeric simulations 
|S] have predicted that prior to the crystallization, in the range of the interaction 
parameter 25 < Vg < 35, the ground state of the system is a strongly correlated 
ferromagnetic Fermi liquid. A paramagnetic Fermi liquid is realized at yet higher 
electron densities when the interactions are relatively weak (r^ ~ 1). The effective 
mass, m, and Lande g factor within the Fermi liquid theory are renormalized due to 
spin exchange effects, with renormalization of the g factor being dominant compared 
to that of the effective mass jS]. Alternatively, near the onset of Wigner crystallization, 
strong increase of the effective mass is expected |3 |S] . 

Recently, there has been a lot of interest in the electron properties of dilute 2D 
systems due to the observation of an unexpected metal-insulator transition (MIT) in 
zero magnetic field, strong metallic temperature dependence of the resistance in these 
systems, and a giant positive magnetoresistance in a magnetic field parallel to the 2D 
plane (for a review, see Ref. The most pronounced effects have been observed in 
high-mobility silicon metal-oxide-semiconductor field-effect transistors (MOSFETs), 
with the low-temperature drop of the resistance reaching an order of magnitude 
and magnetoresistance exceeding five orders of magnitude. Significant progress has 
been recently made in understanding the metallic behaviour of the resistivity and its 
suppression by a magnetic field |1UI the metal-insulator transition, however, still 
lacks adequate theoretical description. 

In this paper, which summarizes results obtained in Refs. |12l 1131 n~Tl 1151 llHj . 
we report measurements of the effective mass, m, and Lande g factor in a wide 
range of electron densities including the immediate vicinity of the MIT. We have 
used two independent methods to determine these parameters. In the first method, 
we have studied low-temperature magnetotransport in a parallel magnetic field. It 
turns out that the magnetic field. Be, required to fully polarize electron spins, is a 
strictly linear function of the electron density: Be oc {rig — n-^) where is some finite 
electron density close to the critical electron density ric for the B = metal-insulator 
transition. Vanishing Be points to a sharply increasing spin susceptibility, x gm, 
and gives evidence in favour of the spontaneous spin polarization at a finite electron 
density. (Similar conclusion about possible spontaneous spin polarization in the dilute 
2D electron system in silicon has been reached in Ref. ^2].) Comparing our data for 
zero- field resistivity with the recent theory JOI) we extract the values of m and g 
separately. It turns out that it is the value of the effective mass that becomes strongly 
(by more than a factor of 3) enhanced with decreasing electron density, while the g 
factor remains nearly constant and close to its value in bulk silicon. In the second 
method, we have determined the effective mass by analyzing temperature dependence 
of the weak-field Shubnikov-de Haas (SdH) oscillations and found good agreement 
with the data obtained by the first method. Furthermore, using tilted magnetic 
fields, we find that the value of the effective mass does not depend on the degree 
of spin polarization, which points to a spin-independent origin of the effective mass 
enhancement. This is in clear contradiction with existing theories [H1[3|S]. 

2. Samples and experimental technique 

Measurements have been made in a rotator-equipped Oxford dilution refrigerator with 
a base temperature of « 30 mK on low-disordered (lOO)-silicon samples similar to 
those previously used in Ref. ^H!- Peak electron mobilities in these samples are close 
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Figure 1. Shubnikov-de Haas oscillations in the Si MOSFET at T Ri 40 mK 
(a) at a relatively high electron density Us = 5.43 X 10^^ cm~^ and (b) at low 
electron density ris = 9.3 X 10^" cm~^. The minima of the resistance at Landau 
level filling factors u = 6 and 10 are shown on an expanded scale in the inset. 



to 3 m^/Vs at 0.1 K. To minimize the contact resistance, which tends to grow very high 
at mK temperatures and low electron density, thin gaps in the gate metallization have 
been introduced; this allows for maintaining high electron density near the contacts 
(about 1.5 X 10"'^^ cm~^) regardless of its value in the main part of the sample. These 
gaps are narrow enough (< 100 nm) for the given gate-oxide thickness to provide a 
smoothly descending electrostatic potential from the high-density part to the low- 
density part. The resistance, Rxx, has been measured by a standard 4-terminal 
technique at a low frequency (0.4 Hz) to minimize the out-of-phase signal. Excitation 
current has been kept low (0.1-0.2 nA) to ensure that measurements are taken in the 
linear regime of response; the power generated in the samples has been maintained 
under 10"^'' W. To verify that the electrons are not overheated, we have studied the 
temperature dependence of the amplitude of the SdH oscillations; the latter has been 
found to follow the theoretical curve down to temperatures less than 50 mK (for more 
on this, see below). 



3. Experimental results 

We start by showing a low-temperature longitudinal magnetoresistance R^x in a 
perpendicular magnetic field B± for a relatively high (Fig. Q^a)) and relatively low 
(Fig. n^b)) electron densities. At the high density, positions of SdH oscillations 
correspond to "cyclotron" filling factors jj^li some of which are marked by arrows. 
Indeed, the energy splittings As = gfiBB± at "spin" filhng factors, v = 2, 6, 
10... = 4i — 2, in high-density Si MOSFETs are known to be much smaller than 
the splittings Ac = hilc — giJ'BB± at "cyclotron" filling factors, v — 8, 12... = 4i, 
disregarding the odd ly valley splitting which is small (here flc — eB±/mc is the 
cyclotron frequency and i = 1, 2, 3...). The behaviour of the sample at a relatively high 
electron density is thus rather ordinary. In contrast, at low electron density (just above 
the metal-insulator transition which in this sample occurs at Us = ric = 8x 10^" cm~^), 
the magnetoresistance looks quite different [20]. The resistance minima are seen only 
at 1^ = 2, 6, and 10 (see the inset); there is also a minimum at v — I (not shown 
in the figure) corresponding to the valley splitting. There are neither dips nor other 
anomalies at magnetic fields corresponding to = 4, 8, or 12 where cyclotron minima 
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Figure 2. Evolution of the Shubnikov-de Haas oscillations with electron density 
in two ranges of filling factors: (a) 3 < u < 9 and (b) 8 < u < 17; T ^ 40 mK. 
The curves are vertically shifted for clarity. 

are expected. 

Figure |21 shows how the resistance minima corresponding to the cyclotron 
spUttings gradually disappear as the electron density is reduced. At the highest 
electron densities (the lower curves), deep resistance minima near even filling factors 
are seen (i/ = 4, 6, and 8 in Fig.|2Ia); — 10, 12, and 16 in Fig.[2Ib)), and a shallow 
minimum is visible at = 14 in Fig. I^^b). As Ug is reduced, the minima aX v — A, 8, 
12, and 16 become less deep, and at the lowest electron densities (the upper curves), 
neither of them is seen any longer, and only minima ai v — 6, 10, and 14 remain. 

Our results thus show that as one approaches the metal-insulator transition, 
the energy gaps at "cyclotron" filling factors become gradually smaller than those 
at "spin" filling factors and eventually vanish. The condition for vanishing Ac is 
gfJ^sB^ = hflc (within the uncertainty associated with the broadening of the energy 
levels), or gm/2me = 1, which is higher by more than a factor of five than the "normal" 
value of this ratio, gm/2me = 0.19. Therefore, the spin susceptibility x oc gm is 
strongly enhanced near the MIT. 

One could attempt to link the observed behaviour to a many-body enhancement of 
spin gaps specific for a perpendicular magnetic field [23 • However, the disappearance 
of the cyclotron splittings in a wide range of magnetic fields would require an enhanced 
^-factor which is independent of magnetic field, in contradiction with Ref. |21| . 
On the other hand, our results are consistent with the suggestion that the 
effective g factor is nearly field-independent and approximately equal to its many-body 
enhanced zero-field value. To probe this conjecture, we have studied the parallel-field 
magnetotransport in a wide range of electron densities. 

Typical curves of the low-temperature magnetoresistance p(i?||) in a parallel 
magnetic field are displayed in Fig. U] Note that the thickness of the 2D electron 
system in Si MOSFETs is small compared to the magnetic length in accessible 
fields, and, therefore, the parallel field couples largely to the electrons' spins while 
the orbital effects are suppressed j2SlE]- The resistivity increases with field until 
it saturates at a constant value above a certain density-dependent magnetic field. 
According to Refs. |^ the saturation of the magnetoresistance indicates the 
onset of a complete spin polarization. In the vicinity of the metal-insulator transition, 
the magnetoresistance is strongly T-dependent down to the lowest achievable 
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Figure 3. Low-temperature magnetoresistance in parallel magnetic fields at 
different electron densities above the critical density for the B = metal-insulator 
transition. The lowest density curve is outside the scaling region as described in 
the text. 

temperatures. As one moves away from the transition, however, the temperature 
dependences saturate at very low temperatures. The data shown below are obtained 
in this low-temperature limit where the magnetoresistance becomes temperature- 
independent. 

In Fig. 21 we show how the normalized magnetoresistance, measured at different 
electron densities, collapses onto a single curve when plotted as a function of B^^/Bc. 
The scaling parameter, Be, has been normalized to correspond to the magnetic field 
at which the magnetoresistance saturates (within the accuracy with which the latter 
can be determined). The observed scaling is remarkably good for B\\/Bc < 0.7 in 
the electron density range 1.08 x 10^^ — 10^^ cm~^, although with increasing n^, 
the scaled experimental data occupy progressively shorter intervals on the resulting 
curve. Both at B\\/Bc > 0.7 and outside the indicated range of electron densities, the 
scaled data start to noticeably deviate from the universal curve. In particular, the 
scaling breaks down when one approaches {us < 1.3 ric) the metal-insulator transition 
which in this sample occurs at zero magnetic field at ric = 8 x 10^*^ cm~^. This is 
not surprising as the magnetoresistance near Uc depends strongly on temperature, as 
discussed above. We note that the observed scaling dependence is described reasonably 
well by the theoretical dependence of p/p{0) on the degree of spin polarization 
^ = gmfiBBw/irh^ns — B\\/Bc predicted by the recent theory |26| . 

In Fig. 13 Be is plotted vs. Us- With high accuracy. Be is proportional to the 
deviation of the electron density from its critical value, i.e., to (n^ — ric), over a wide 
range of electron densities. In other words, the field, at which the magnetoresistance 
saturates, tends to vanish at tIc (see also Refs. ^| EZ|)- We emphasize that 
our procedure provides high accuracy for determining the behaviour of the field of 
saturation with electron density, i.e., the functional form of Bc{ns), even though 
the absolute value of Be is determined not so accurately. Note that at Us above 
2.4 X 10^^ cm~^, the saturation of the resistance is not reached in our magnetic field 
range; still, the high precision of the collapse of the high-density experimental curves 
onto the same scaling curve as the low-density data allows us to draw conclusions 
about the validity of the obtained law Bc(ns) over a much wider range of electron 
densities. 
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Figure 4. Scaled curves of the normalized magnetoresistance at different Us vs. 
B^^/Bc. The electron densities arc indicated in units of 10^^ cm~^. Also shown 
by a dashed line is the normalized magnetoresistance calculated in Ref. 1261 . 
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Figure 5. Dependence of the field Be on electron density. The dashed line 
is a linear fit which extrapolates to the critical electron density for the B = 
metal-insulator transition. Comparison of l/A*(ns) (see Eq. 0) and Bc{ns) 
(open and solid circles, correspondingly) is shown in the inset. The dashed lines 
are linear fits which extrapolate to the critical electron density for the metal- 
insulator transition. 



The observed tendency of Be to vanish at a finite electron density is consistent 
with the strong increase of the spin susceptibihty x oc gm [12i 1151 I28| and gives 
evidence in favour of a spontaneous spin polarization at n^^ Pi ric- In principle, either 
g 01 m (or both) may be responsible for the strong increase of the spin susceptibility. 
As has already been mentioned, within the Fermi liquid theory, both the effective 
mass and g factor are renormalized due to spin exchange effects, with renormalization 
of the g factor being dominant compared to that of the effective mass. In contrast, 
the dominant increase of the effective mass follows from an alternative description of 
the strongly- interacting electron system beyond the Fermi liquid approach [3 El ■ To 
separate g and m, we have measured the temperature-dependent conductivity in zero 
magnetic field and analyzed the data in the spirit of recent theory (10| . According to 
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Figure 6. The temperature dependence of the normaUzed conductivity at 
different electron densities (indicated in units of lO^^cm"'^) above the critical 
electron density for the metal-insulator transition. The dashed lines are fits of 
the linear interval of the dependence. 



this theory, ct is a hnear function of temperature: 
a{T) 



= 1 - A*kBT, 



(1) 



where the slope, A*, is determined by the interaction-related parameters: the Fermi 
liquid constants, Fq and : 



A* 



(1 + aF^)gm 



(2) 



The factor a is equal to 8 in our case [221 • This theoretical relation allows us to 
determine the many-body enhanced g factor and mass m separately using the data 
for the slope A* and the product gm. 

Typical dependences of the normalized conductivity on temperature, a{T)/aQ, 
are displayed in Fig. at different electron densities above ric, the value cto, which 
has been used to normalize cr, was obtained by extrapolating the linear interval of the 
cr(r) dependence to T = 0. As long as the deviation \a/ao — l\ is sufficiently small, the 
conductivity a increases linearly with decreasing T in agreement with Eq. until 
it saturates at the lowest temperatures. As seen from the figure, the linear interval of 
the dependence is wide enough to make a reliable fit. 

The Us dependence of the inverse slope 1/A*, extracted from the a{T) data, is 
shown in the inset to Fig. [S] by open circles. Over a wide range of electron densities, 
the values 1/A* and ^bBc turn out to be close to each other. The low-density data for 
1/A* are approximated well by a linear dependence which extrapolates to the critical 
electron density ric in a similar way to the behaviour of the polarization field Be. 

In Fig. 13 we show the so-determined values 5/50 and m/rrib as a function of the 
electron density (here 50 = 2 is the g factor in bulk silicon, mb is the band mass equal 
to 0.19r7ie, and TOg is the free electron mass). In the high region (relatively weak 
interactions), the enhancement of both g and m is relatively small, both values slightly 
increasing with decreasing electron density in agreement with earlier data Also, 
the renormalization of the g factor is dominant compared to that of the effective mass, 
which is consistent with theoretical studies 0. 
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Figure 7. The effective mass (circles) and g factor (squares), determined from 
the analysis of the parallel field magnetoresistance and temperature-dependent 
conductivity, vs. electron density. The dashed lines are guides to the eye. The 
inset compares the so-obtained effective mass (dotted line) with the one extracted 
from the analysis of SdH oscillations (circles). 



In contrast, the renormalization at low (near the critical region), where ^ 1, 
is much more striking. As the electron density is decreased, the renormalization of 
the effective mass overshoots abruptly while that of the g factor remains relatively 
small, g ~ go, without tending to increase. Hence, the current analysis indicates that 
it is the effective mass, rather than the g factor, that is responsible for the drastically 
enhanced gm value near the metal-insulator transition. 

Since the procedure for extracting g and m described above relies on theoretically 
calculated functional form for the slope A* we have performed independent 

measurements of the effective mass based on the temperature analysis of the 
amplitude. A, of the weak-field (sinusoidal) SdH oscillations. A typical temperature 
dependence of A for the normalized resistance, Rxx/Ro (where Rq is the average 
resistance), is displayed in Fig.|S| To determine the effective mass, we use the method 
of Ref. |31| extending it to much lower electron densities and temperatures. We fit 
the data for A{T) using the formula 

where Aq — 4exp(— 27r^fcBTD/?iilc) and To is the Dingle temperature. As the latter 
is related to the level width through the expression To — h/2TTkBT (where r is the 
elastic scattering time) damping of the SdH oscillations with temperature may 
be influenced by temperature-dependent r. We have verified that in the studied low- 
temperature limit for electron densities down to w 1 x 10^^ cm~^, possible corrections 
to the mass value caused by the temperature dependence of r (and hence To) are 
within our experimental uncertainty which is estimated at about 10%. Note that 
the amplitude of the SdH oscillations follows the calculated curve down to the lowest 
achieved temperatures, which confirms that the electrons were in a good thermal 
contact with the bath and were not overheated. The fact that the experimental 
dependence A{T) follows the theoretical curve justify applicability of Eq. Q to this 
strongly-interacting electron system. 
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Figure 8. Amplitude of the weak-field SdH oscillations vs. temperature at 
Us = 1.17 X 10^^ cm~^ for oscillation numbers u = hcns/eB±^ = 10 (dots) and 
u = lA (squares). The value of T for the v = 10 data is divided by the factor of 
1.4. The solid line is a fit using Eq. ISl . 
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Figure 9. The effective mass vs. the degree of spin polarization for the 
following electron densities in units of 10^^ cm~^: 1.32 (dots), 1.47 (squares), 
2.07 (diamonds), and 2.67 (triangles). The dashed lines are guides to the eye. 

The so-determined effective mass is shown in the inset to Fig. [7\ In quantitative 
agreement with the results obtained by the alternative method described above (the 
dotted line), the effective mass sharply increases with decreasing n^. The agreement 
between the results obtained by two independent methods adds confidence in our 
results and conclusions. Our data are also consistent with the data for spin and 
cyclotron gaps obtained by magnetocapacitance spectroscopy 22i- 

A strong enhancement of m at low electron densities may originate from spin 
effects iniEIISI- With the aim of probing a possible contribution from the spin effects, 
we have introduced a parallel magnetic field component to align the electrons' spins. In 
Fig.O we show the behaviour of the effective mass with the degree of spin polarization, 
p = {B^ -|- Bp^/^/Bc. As seen from the figure, within our accuracy, the effective mass 
m does not depend on p. Therefore, the m{ns) dependence is robust, the origin of the 
mass enhancement has no relation to the electrons' spins and exchange effects 
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4. Discussion 

Under the conditions of our experiments, the interaction parameter, rg, is larger by a 
factor of 2m/mb than the Wigner-Seitz radius and reaches approximately 50, which is 
above the theoretical estimate for the onset of Wigner crystallization. As has already 
been mentioned, two approaches to calculate the renormalization of m and g have been 
formulated. The first one exploits the Fermi liquid model extending it to relatively 
large Tg. Its main outcome is that the renormalization of g is large compared to that of 
m In the limiting case of high r^, one may expect a divergence of the g factor that 
corresponds to the Stoner instability. These predictions are in obvious contradiction 
to our data: (i) the behaviour of the 2D dilute system in the regime of the strongly 
enhanced susceptibility — close to the onset of spontaneous spin polarization and 
Wigner crystallization — is governed by the effective mass, rather than the g factor, 
through the interaction parameter r^; and (ii) the insensitivity of the effective mass 
to spin effects also cannot be accounted for. 

The other theoretical approach either employs analogy between a strongly 
interacting 2D electron system and He'^ |7] or applies Gutzwiller's variational method 

to Si MOSFETs 8 . It predicts that near the crystallization point, the 
renormalization of m is dominant compared to that of g and that the effective mass 
tends to diverge at the transition. Although the sharp increase of the mass is in 
agreement with our findings, it is the expected dependence of m on the degree of 
spin polarization that is not confirmed by our data: the model of Ref. [7] predicts 
that the effective mass should increase with increasing spin polarization, whereas the 
prediction of the other model ^ is the opposite. 

Thus, the existing theories fail to explain our finding that in a dilute 2D electron 
system the effective mass is strongly enhanced and does not depend on the degree of 
spin polarization. The fact that the spin exchange is not responsible for the observed 
mass enhancement reduces the chances for the occurrence of the ferromagnetic Fermi 
liquid prior to the Wigner crystallization. However, should the spin exchange be small, 
the spin effects may still come into play closer to the onset of Wigner crystallization 
where the Fermi energy may continue dropping as caused by mass enhancement. 

In summary, we have found that in very dilute two-dimensional electron systems 
in silicon, the effective mass sharply increases with decreasing electron density, while 
the g factor remains nearly constant and close to its value in bulk silicon. The 
enhanced effective mass does not depend on the degree of the spin polarization and, 
therefore, its increase is not related to spin exchange effects, in contradiction with 
existing theories. The corresponding strong rise of the spin susceptibility may be a 
precursor of a spontaneous spin polarization; unlike in the Stoner scenario, the latter 
originates from the enhancement of the effective mass rather than the increase of the 
g factor. Our results show that the dilute 2D electron system in silicon behaves well 
beyond the weakly interacting Fermi liquid. 
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